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Influence of cavity design 
preparation on stress values in 
maxillary premolar: 
a finite element analysis 



Aim To analyze the influence of cavity design preparation 
on stress values in three-dimensional (3D) solid model of 
maxillary premolar restored with resin composite. 

Methods 3D solid model of maxillary second premolar 
was designed using computed-tomography (CT) data. 
Based on a factorial experiment, 9 different mesio-occlus- 
al-distal (MOD) cavity designs were simulated, with three 
cavity wall thicknesses (1.5 mm, 2.25 mm, 3.0 mm), and 
three cusp reduction procedures (without cusp reduction, 
2.0 mm palatal cusp reduction, 2.0 mm palatal and buc- 
cal cusp reduction). All MOD cavities were simulated with 
direct resin composite restoration (Gradia Direct Posterior, 
GC, Japan). Finite element analysis (FEA) was used to calcu- 
late von Mises stress values. 

Results The von Mises stresses in enamel, dentin, and 
resin composite were 79.3-233.6 MPa, 26.0-32.9 MPa, and 
1 80.2-252.2 MPa, respectively. Considering the influence of 
cavity design parameters, cuspal reduction (92.97%) and 
cavity wall thickness (3.06%) significantly (P<0.05) deter- 
mined the magnitude of stress values in enamel. The influ- 
ence of cavity design parameters on stress values in dentin 
and resin composite was not significant. When stresses for 
enamel, dentine, and resin composite were considered all 
together, palatal cusp coverage was revealed as an optimal 
option. Cavity wall thickness did not show a significant ef- 
fect on stress values. 

Conclusion Based on numerical simulations, a palatal cusp 
reduction could be suggested for revealing lower stress 
values in dental tissues and restorative material. This type 
of cavity design should contribute to better biomechanical 
behavior of tooth-restoration complex, consequently pro- 
viding the long-lasting clinical results. 
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in the recent years, there has been an increasing interest 
in the research of biomechanical aspects of biomaterials 
and human tissues (1-3). Although studies conducted in 
vivo and in vitro have provided some of the answers in 
this field, dental and medical research is usually costly, and 
may be ethically questionable and time-consuming (4,5). 
Because of this, the use of numerical models and in vitro 
simulations became a valuable tool for saving time and 
money associated with laboratory and clinical research 
(6). Previous studies have reported different techniques 
for generating three-dimensional (3D) solid models of the 
teeth (7-9). Nowadays, technological development brings 
new possibilities for efficient generation of sophisticated 
3D solid models. For example, using specialized software, 
these models can be generated based on computed-to- 
mography (CT) scan data (5,10-12). In addition, the appli- 
cation of finite element analysis (FEA) allows calculation of 
stress and strain within tooth structure and biomaterials, 
which can hardly be measured in vivo (1 3). 

Cavity design preparation has a great impact on stress val- 
ues and fracture resistance of a tooth (7,1 4-16). It is a factor 
of a paramount importance, especially in cases of restoring 
maxillary premolars with extensive mesio-occlusal-distal 
(IVIOD) cavities (14,1 7,18). From the biomechanical point of 
view, different opinions have been reported on the most 
appropriate restorative procedure in such cases. Kuijs et al 
have found that ceramic, indirect resin composite and di- 
rect resin composite provide comparablefatigue resistance 
in a cusp replacing restorations (19). These findings were 
supported by clinical trials performed by van Dijken and 
Hickel et al (20,21). On the other hand, Soares et al found 
that MOD cavities restored with resin composite placed 
with direct technique attained better biomechanical per- 
formance than those restored with laboratory processed 
resin and ceramic restorations (15,16). Another study also 
confirmed that in comparison with ceramic restorations, 
resin composite restoration had higher fatigue resistance 
(17). As opposed to preparation for direct restoration, cav- 
ity preparation for indirect restorations requires removal of 
additional amount of tooth structure (22). The situation is 
the same with the cavity preparation for the amalgam (23). 
Since the quantity of the tooth structure removed while 
doing cavity preparation affects the biomechanical char- 
acteristics of the restored tooth, the use of adhesive direct 
restorations should be recommended for reinforcing the 
remaining dental structure (2,16). 

When planning the design of MOD cavity preparation, cav- 
ity wall thickness and cusp reduction should be carefully 



considered. Usually, cusp reduction is recommended when 
cavity isthmus width is 2/3 of intercuspal width (7,24). Al- 
though this promotes more dental tissue reduction (25), it 
was shown that the reduction of cuspal height by 2.0 mm 
increases fracture resistance of a premolar when restored 
with direct resin composite (26,27). On the other hand, 
cavity wall thickness is not well defined. Macpherson et al 
found that 2.25 mm wall thickness is critical for restoring 
fracture resistance of tooth with MOD cavity (28), but an- 
other study, which investigated 1 .0-3.0 mm wall thickness- 
es, reported that the thickness of remaining cavity walls 
was not relevant to fracture resistance (27). 

The aim of this study was to investigate the effect of cav- 
ity design preparation on stress values in remaining tooth 
structures restored with resin composite. The null hypoth- 
esis was that stress values were not significantly influenced 
by the cavity wall thickness and cusp reduction. 

MATERIAL AND METHODS 

The study was conducted at the Clinic of Dentistry of 
Vojvodina, Novi Sad, Serbia, and at the Institute of Physics, 
Belgrade, Serbia, from December 201 1 to March 2012. A 
human maxillary second premolar, extracted for orthodon- 
tic reasons, at the Clinic of Dentistry of Vojvodina was used 
for the study. Immediately after the extraction, the tooth 
was cleaned of soft tissue remnants and used for 3D solid 
model generation. The selected tooth was intact, without 
caries, fractures, and morphological abnormalities. 

3D solid model generation 

The extracted tooth was scanned using multilayer CT scan- 
ner (SOMATOM Sensation 64 Cardiac, Siemens, Forchheim, 
Germany). A total of 1 1 0 slices were made along x-axis, 88 
along y-axis, and 47 along z axis. For the solid model gener- 
ation, slices along z axis were used with 0.5 mm resolution. 
The selected slices were imported to AMIRA software (Vis- 
age Imaging Inc, San Diego, CA, USA) for automatic tooth 
structures (enamel, dentin, pulp) segmentation. The seg- 
mentation was based on image density threshold of differ- 
ent gray scale intensities corresponding to various degrees 
of mineralization. Obtained contours were then imported 
into SolidWorks 2011 software (Dassault Systcmes Solid- 
Works Corp, Waltham, MA, USA), and 3D solid model of the 
intact maxillary second premolar was generated by using 
a lofting technique. Additionally, based on the outer ge- 
ometry of the model, periodontal ligament and alveo- 
lar bone (cortical and cancellous) were created. 
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Cavity preparation design 

In the 3D solid model of intact maxillary premolar, differ- 
ent MOD cavities were designed. All MOD cavities had pul- 
pal and axial walls with at least 1.0 mm dentin thickness 
over the pulp (29), while the gingival walls were located 
1 .0 mm above cemento-enamel junction (CEJ). Solid mod- 
els of MOD cavities were created with three different wall 
thicknesses (3.0 mm, 2.25 mm, and 1 .5 mm), and three dif- 
ferent cusp reduction procedures (without cusp reduc- 
tion, 2.0 mm palatal cusp reduction, and 2.0 mm palatal 
and buccal cusp reduction). In total, there were nine 3D 
solid tooth models (Figure 1). Restorations of all MOD cavi- 




FIGURE 1 . 3D solid tooth models with different mesio-occlusal- 
dlstal (MOD) cavities. (A) Model 1 : 3D solid tooth model with 
MOD cavity with 3.0 mm wall thickness, without cusp reduc- 
tion; (B) Model 2: 3D solid tooth model with MOD cavity with 
2.25 mm wall thickness, without cusp reduction; (C) Model 
3: 3D solid tooth model with MOD cavity with 1 .5 mm wall 
thickness, without cusp reduction; (D) Model 4: 3D solid tooth 
model with MOD cavity with 3.0 mm wall thickness and 2.0 
mm palatal cusp reduction; (E) Model 5: 3D solid tooth model 
with MOD cavity with 2.25 mm wall thickness and 2.0 mm 
palatal cusp reduction; (F) Model 6: 3D solid tooth model with 
MOD cavity with 1 .5 mm wall thickness and 2.0 mm palatal 
cusp reduction; (G) Model 7: 3D solid tooth model with MOD 
cavity with 3.0 mm wall thickness and 2.0 mm palatal and 
buccal cusp reduction; (H) Model 8: 3D solid tooth model with 
MOD cavity with 2.25 mm wall thickness and 2.0 mm palatal 
and buccal cusp reduction; (I) Model 9; 3D solid tooth model 
with MOD cavity with 1 .5 mm wall thickness and 2.0 mm 
palatal and buccal cusp reduction. 



ties were simulated as direct resin composite restorations 
(Gradia Direct Posterior, GC, Tokyo, Japan). 

Finite element analysis 

Nine 3D solid tooth models with MOD cavities restored 
with resin composite were used to simulate nine clinically 
different cavity designs. All 3D solid models were derived 
from the 3D solid model of the intact maxillary second pre- 
molar 3D solid models were meshed with parabolic tet- 
rahedral elements. The parabolic tetrahedral element is 
defined by four corner nodes, six mid-side nodes, and six 
edges. These elements were used because they represent 
curved boundaries more accurately and provide better 
mathematical approximations. The number of elements 
and nodes varied according to the model (142,407-1 75 727 
elements and 223113-268918 nodes). Convergence test 
was used to verify that our numerical model reached con- 
verged results and that no further mesh refinement was 
necessary. The exterior nodes on all surfaces of the cortical 
bone were restrained in all directions as the boundary con- 
ditions for all models. All materials were assumed to have 
linear, elastic, and isotropic properties (7,8,10), represented 
by the Young's modulus of elasticity and the Poisson's ratio 
(Table 1) (7,14,30). 

To simulate masticatory forces in the maximum intercuspa- 
tion position, the static axial load with the resulting force of 
200 N was applied on the occlusal surface of a tooth at three 
points (palatal cusp tip and both marginal ridges) (31,32). 
Using structural FEA, maximum von Misses stress values in 
the enamel, dentin, and resin composite were calculated. 

Statistical methods 

To investigate the impact of cavity design parameters on 
stress values in maxillary premolar, a full three-level factori- 
al design based on a quadratic model was used with a total 
of nine (3^) experiments (33). Two design factors were con- 



TABLE 1 . Material properties assigned to dental tissues and 
restorative material 





Young's 


Poisson's 




Material 


modulus (MPa) 


ratio 


References 


Enamel 


84100 


0.20 


(714) 


Dentin 


18600 


0.31 


(714) 


Pulp 


2 


0.45 


(714) 


Periodontal ligament 


70 


0.45 


(714) 


Cortical bone 


15000 


0.30 


(714) 


Cancellous bone 


1500 


0.30 


(714) 


Resin composite 


6700 


0.22 


(30) 
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sidered: cavity wall thickness and cuspal reduction. Each 
design factor was assigned three levels (Table 2). 

In order to establish the relative importance of the inves- 
tigated factors and their interactions, analysis of variance 
(ANOVA) was conducted for enamel, dentin, and resin 
composite stress values, respectively (7). Factorial analysis 
was conducted in Statistica vlO software (Statsoft Inc, Tul- 
sa, OK, USA). 



(30.80%) affected the magnitude of stress values (Table 
5). However, no significant association (P<0.05) was es- 
tablished. The main effects plot for dentin stress again re- 
vealed that the increase of cavity wall thickness resulted in 
lower stresses, while in the case of cuspal reduction, pala- 
tal cusp reduction resulted in lower stress values (Figure 4), 
The interaction plot (Figure 5) indicated that, at maximum 
cavity wall thickness, the application of palatal and buccal 
cusp reduction resulted in the lowest stress values. 



RESULTS 



Main Effects Plot for Enamel stress 



The von Mises stresses in enamel, dentin, and resin com- 
posite, obtained by FEA, were 79.3-233.6 IVlPa, 26.0-32.9 
IVlPa, and 180.2-252.2 MPa, respectively (Table 3). 

Considering the influence of cavity design parameters on 
enamel stress, the ANOVA showed that cuspal reduction 
(92.97%) and cavity wall thickness (3.06%) significantly 
(P<0.05) determined the magnitude of stress values (Ta- 
ble 4). The main effects plot for enamel stress (Figure 2) 
revealed that the increased cavity wall thickness moder- 
ately reduced enamel stress, while the presence of palatal 
or palatal and buccal cusp reduction contributed to lower 
stress values, compared to the case with no cuspal reduc- 
tion. Interaction plot (Figure 3) showed that at maximum 
wall thickness value, application of palatal and buccal cusp 
reduction resulted in minimal stress values. 

For the stress values of dentin, the results also showed 
that cavity wall thickness (34.92%) and cuspal reduction 

TABLE 2. Factorial design - investigated factors and assigned 
levels 

Level 




1,5 3.0 none palatal pal_buccal 

Cavity wall thickness [mm] Cuspai reduction 

FIGURE 2. Main effects plot for enamel stress. 

Interaction Plot for Enamel stress 



240 " Cuspal reduction=none 



200 



i: 160 - 



E 

" 120 - 



Investigated factor 


1 


2 


3 


Cavity wall thickness (mm) 


1.50 


2.25 


3.00 


Cuspal reduction (mm) 


None 


Palatal 


Palatal and 






(2.00) 


buccal (2.00) 



Cuspal 
80 <— 



Cuspal reduct.=palatal 
reiduct.=pal_buccal 



' Cuspal reiduction=none 



Cuspal reduct.=pal_buccaL 
Cuspal reducl.=palatal 



1.5 3.0 
Cavity wall thickness [mm] 



FIGURE 3. Interaction plot for enamel stress. 
TABLE 3. Factorial experiments in randomized order, with respective output stresses obtained by finite element analysis 



Run 


Cavity wall thickness (mm) 


Cuspal reduction 


Enamel stress (MPa) 


Dentin stress (MPa) 


Resin composite stress (MPa) 


1 


2.25 


palatal 


79.3 


26,0 


180.8 


2 


1.50 


palatal and buccal 


83.4 


32.8 


199.7 


3 


3.00 


palatal and buccal 


82.7 


27.2 


196.6 


4 


2.25 


none 


196.6 


32.9 


188.9 


5 


2.25 


palatal and buccal 


878 


26,9 


182.5 


6 


T50 


none 


233.6 


32.1 


252.2 


7 


T50 


palatal 


94.2 


28,7 


180.2 


8 


3.00 


palatal 


83.2 


27.5 


195.2 


9 


3.00 


none 


172.5 


27,7 


210.3 
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Regarding the influence of cavity design on restorative 
material, tine analysis of stress values in resin compos- 
ite revealed the dominant effect of cuspal reduction 
(40.87%), followed by interaction of cavity wall thickness 
and cuspal reduction (20.94%), and cavity wall thickness 
(3.72%) (Table 6). However, as in the case of dentin stress. 

Main Effects Plot for Dentin stress 



32 : 




26 b 

1,5 3.0 none palatal pal_buccal 

Cavity wall thickness [mm] Cuspal reduction 

FIGURE 4. Main effects plot for dentin stress. 



no significant association (P<0.05) was established. The 
main effects plot showed that the minimum stress was 
obtained for medium cavity wall thickness, ie, the palatal 
cusp reduction (Figure 6). From the interaction plot (Fig- 
ure 7) it is also evident that at medium cavity wall thick- 
ness, the most favorable stress was obtained for palatal 
cusp reduction. 

Stress distribution patterns were similar for all models (Fig- 
ure 8). The highest stress values occurred at loading points 
(palatal cusp tip and both marginal ridges) and at cervical 
area of the palatal surface. 

DISCUSSION 

The present study indicated that cavity design preparation 
generally affected von Mises stress values in premolar re- 
stored with direct resin composite. While the cusp reduc- 
tion decreased stress values and provided more favorable 
stress distribution, cavity wall thickness showed no signifi- 
cant influence. 



TABLE 4. ANOVA statistical results of von Mises stresses for enannel 



Source 


Degrees of 
freedonn 


Sum of 
squares 


Mean sum 
of squares 


Total sum of 
squares (%) 


P 


Cavity wall thickness 


1 


883.31 


883.31 


3.06 


0.0395 


Cuspal reduction 


2 


26820.60 


13410.30 


92.97 


0.0028 


Cavity wall thickness x cavity wall thickness 


1 


2738 


2738 


0.10 


0.4811 


Cavity wall thickness x cuspal reduction 


2 


1044.04 


522.02 


3.61 


0.0664 


Residual 


2 


74,30 


3715 


0.26 




Total 


8 


28849.63 




100.00 




TABLE 5. ANOVA statistical results of von Mises stresses for dentin 


Source 


Degrees of 
freedom 


Sum of 
squares 


Mean sum 
of squares 


Total sum of 
squares (%) 


P 


Cavity wall thickness 


1 


20.91 


20,91 


34.92 


0.2291 


Cuspal reduction 


2 


18.44 


9.22 


30.80 


04362 


Cavity wall thickness x cavity wall thickness 


1 


1.08 


1.07 


1.80 


0.7353 


Cavity wall thickness x cuspal reduction 


2 


5.17 


2.59 


8.64 


0.7339 


Residual 


2 


14.27 


714 


23.84 




Total 


8 


59.87 




100.00 





TABLE 6. ANOVA statistical results of von Mises stresses for resin composite 



Source 


Degrees of 
freedom 


Sum of 
squares 


Mean sum 
of squares 


Total sum of 
squares (%) 


P 


Cavity wall thickness 


1 


150.00 


150.00 


3.72 


0.5019 


Cuspal reduction 


2 


164940 


824.70 


40.87 


0.2161 


Cavity wall thickness x cavity wall thickness 


1 


936.00 


936.00 


23.20 


0.1796 


Cavity wall thickness x cuspal reduction 


2 


845.11 


422,56 


20.94 


0.3498 


Residual 


2 


454.70 




11.27 




Total 


8 


4035.21 




100.00 
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Interaction Plot for Dentin stress 






Interaction Plot for Resin composite stress 
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FIGURE 5. Interaction plot for dentin stress. 



FIGURE 7. Interaction plot for resin composite stress. 



iVIain Effects Plot for Resin composite stress 
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Q. 
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in 
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Q. 

E 
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c 
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none palatal pal_buccal 
Cuspal reduction 



FIGURE 6. Main effects plot for resin composite stress. 

Preservation of sound tooth structures is the primary goal of 
modern restorative dentistry. However, from biomechanical 
point of view, protection of remaining tooth structures from 
unfavorable mechanical responses should be considered a 
priority, even if it requires the removal of additional dental 
tissue (7,17). Due to their unfavorable anatomy, maxillary 
premolars with extensive IVIOD cavities are at great risk of 
fracturing if restored without regarding protective principles 
(1,22). It is well known that palatal cusp of maxillary premo- 
lars fractures more frequently than buccal cusp (7,34). Also, 
Soares et al experimentally confirmed that palatal cusp was 
involved in more severe type of fracture than buccal cusp 
(15). In order to prevent this, several studies recommended 
cusp reduction from at least 1.5-3.5 mm (7,17,26,27). In this 
way, higher fracture resistance of restored tooth can be ac- 
complished (17,26,27). Even if fracture occurs, it is going to 
be less severe than without cusp reduction. Some studies 
support the coverage of palatal cusp only (1,27,35), while 
others assert coverage of both cusps (1 7,26). 




stress (MPa) 



FIGURE 8. Von iVlises stress distribution results (A) Model 1 : 
3D solid tooth model with mesio-occlusal-distal (MOD) cavity 
with 3.0 mm wall thickness, without cusp reduction; (B) Model 
2: 3D solid tooth model with MOD cavity with 2.25 mm wall 
thickness, without cusp reduction; (C) Model 3: 3D solid tooth 
model with MOD cavity with 1 .5 mm wall thickness, without 
cusp reduction; (D) Model 4: 3D solid tooth model with MOD 
cavity with 3.0 mm wall thickness and 2.0 mm palatal cusp 
reduction; (E) Model 5: 3D solid tooth model with MOD cavity 
with 2.25 mm wall thickness and 2.0 mm palatal cusp reduc- 
tion; (F) Model 6: 3D solid tooth model with MOD cavity with 
1 .5 mm wall thickness and 2.0 mm palatal cusp reduction; (G) 
Model 7: 3D solid tooth model with MOD cavity with 3.0 mm 
wall thickness and 2.0 mm palatal and buccal cusp reduction; 
(H) Model 8: 3D solid tooth model with MOD cavity with 2.25 
mm wall thickness and 2.0 mm palatal and buccal cusp reduc- 
tion; (I) Model 9: 3D solid tooth model with MOD cavity with 
1 .5 mm wall thickness and 2.0 mm palatal and buccal cusp 
reduction. 



www.cmj.hr 



5?4 



MEDICAL RESEARCH IN BIOPHYSICS 



Croat Med J. 2D12;53:568-?B 



In the present study, the cusp reduction in 3D solid mod- 
els decreased stress values, but significant relationship was 
found only for enamel. The obtained stress values were 
much lower than those calculated for the enamel of intact 
tooth model (155 MPa). Although these results indicated 
similar effect of palatal, and both palatal and buccal cusp 
height reduction, when von Mises stresses for the enamel, 
dentine, and resin composite were considered in conjunc- 
tion, palatal cusp coverage was proven to be the optimal 
option. Stress distribution patterns confirmed this finding. 
Although location of highest stress values was similar for 
all models, stress concentration areas at palatal cusp tip 
for models without cusp reduction were wider than for 
the models with cusp coverage, especially if tooth-res- 
toration interface was in contact with applied force. Fur- 
ther, in models with cusp reduction, stress concentration 
areas were relocated from the remaining tooth structures 
to the restoration. This should preserve remaining tooth 
structures from unfavorable mechanical responses, which 
could make tooth non-restorable (7). Also, no differences 
were found in stress distribution at buccal cusp. This, again, 
implies that only palatal cusp should be reduced in wide 
MOD cavities in maxillary premolars. 

Regarding the cavity wall thickness, the results showed 
that this parameter of cavity design was not as relevant as 
the cusp reduction, which is in accordance with the find- 
ings of other authors (1 4,27). 

In the present study, direct resin composite restoration was 
simulated for all types of cavity preparations. The use of this 
material provides 83.3% fracture resistance of sound tooth 
(15), and stress distribution similar to that of a sound tooth 
(2,16). Further, cavity preparation for direct resin compos- 
ite restoration requires less tooth structure removal than 
for the amalgam or indirect ceramic restoration, which is in 
accordance with the main principle of modern restorative 
dentistry. Also, resin composites have good clinical survival 
rate (90% after two years and 55.1 %-89.7% after ten years), 
even when they are used for restoration of extensive cavi- 
ties in posterior teeth (36,37). 

The 3D solid tooth model was generated using CT scan- 
based data. It presents a modern approach to achieving a 
highly detailed 3Dfinite element model of a tooth (5,10,38), 
and has an important role in investigations of different 
clinical situations in dentistry (8,12). Numerical modeling 
and simulation is useful for obtaining information about 
mechanical behavior of sound and restored tooth (39), 
and is able to demonstrate the otherwise inacces- 



sible stress distribution within the tooth-restoration com- 
plex (1 1 ). It also saves time and costs related to experimen- 
tal studies and clinical trials (5). Furthermore, this method 
allows an infinite number of variables to be studied (6). On 
the other hand, it is obviously quite impossible to include 
all of the factors from the oral environment in a computer 
simulation (40). Moreover, it is known that factors such as 
material properties and loading conditions significantly in- 
fluence the FEA results (14,40). Thus, numerical prototyp- 
ing seems to be a valid method for bringing a new idea 
from concept to clinical application (11). 

In conclusion, based on the numerical simulations and 
analysis applied in this study, a palatal cusp reduction 
could be suggested in order to reveal a lower stress values 
in dental tissues and restorative material. With the respect 
to the sound tooth structures, this type of cavity design 
should contribute to better biomechanical behavior of a 
tooth-restoration complex, consequently providing the 
long-lasting clinical results. 
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